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^ I Ultra-luminous infrared galaxies are among the most luminous objects in the local uni- 

^1 verse and are thought to be powered by intense star formation^ •'^ . It has been shown that 

in these objects the rotational spectral lines of molecular hydrogen observed at mid-infrared 

wavelengths are not affected by dust obscuration^ , leaving unresolved the source of exci- 

CN I tation of this emission. Here I report an analysis of archival Spitzer Space Telescope data 

on ultra-luminous infrared galaxies and demonstrate that star formation regions are buried 

w I inside optically thick clouds of gas and dust, so that dust obscuration affects star-formation 

indicators but not molecular hydrogen. I thereby establish that the emission of H2 is not 

co-spatial with the buried starburst activity and originates outside the obscured regions. 

_. _ This is rather surprising in light of the standard view that H2 emission is directly associated 

^ I with star-formation acti vity^'^'^ . Instead, I propose that H2 emission in these objects traces 

j^ I shocks in the surrounding material, which are in turn excited by interactions with nearby 

galaxies, and that powerful large-scale shocks cooling by means of H2 emission may be much 

more common than previously thought. In the early universe, a boost in H2 emission by 

(-\| ■ this process may speed up the cooling of matter as it collapsed to form the first stars and 

cn ■ 

00 

The analysis is based on Spitzer Space Telescope Infrared (IR) Spectrograph^ mid- 

O ■ IR (5— 35/ini) spectra of 48 ultra-luminous infrared galaxies^ (ULIRGs) publicly available 

^-^ ■ through Spitzer Archive. The objects were selected from the IRAS 1 Jy sample to include 

only those ULIRGs that did not have signatures of powerful supermassive black hole activity 

i^ in their optical spectra^. For these objects I measure the strength of the lO/xm silicate 

^ • absorption feature (S'[9.7/im]), the strength of 6/im water ice absorption {S[6.0fiin]), the 

fluxes of II2 emission lines, the fluxes of polycyclic aromatic hydrocarbon (PAH) emission 

features - which are used as star-formation indicators^ - and the bolometric fluxes F(IR) of 

the galaxies (as described in Supplementary Information). In addition, I use other samples 

of nearby ULIRGs for which some of these measurements are available in the literature'^'^'^'^^. 

As a comparison data set, I use the publicly available spectroscopic data from the Spitzer 

Infrared Nearby Galaxy Survey (SINGS^^), which I analyze using the same methods as 

ULIRG spectra. Nuclei of SINGS galaxies that show activity associated with accretion onto 

a supermassive black hole have been excluded from the comparison sample. Therefore, star 

formation is the likely dominant energy source for both ULIRGs and comparison galaxies 

selected for this Letter. 



galaxies and would make these first structures more readily observable^ . 
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Even though the starbursts supply most of the energy output of ULIRGs, which is in 
excess of 10^^ solar luminosities, only a small fraction of each ULIRG's bolometric luminosity 
is observed as UV and optical continuum of young stars^, as most of this emission is absorbed 
by dust and re-emitted thermally at IR wavelengths. If the observed PAH and H2 emission 
is powered by the embedded starbursts, then these features should also be extincted by 
the intervening matter. Because the amount of extinction depends on wavelength, different 
features will be affected by different amounts (Figured]); the ones that fall within the strong 
opacity peak from astronomical silicates at 9.7/im should be affected most strongly. By 
taking ratios of emission line fluxes within and outside the absorption feature one can test 
where the emitting region is located relative to the source of opacity. Since ULIRGs are 
optically thick even at mid-IR wavelengths, with a median silicate absorption strength of 
1.6, the effect of opacity is expected to be strong. 

Indeed, the observed PAH ratios depend on the strength of absorption (Figure [2]), con- 
firming that the PAH-emitting regions are located behind silicates and water ices. This ob- 
servation provides a quantitative measure of absorption correction that needs to be applied 
to PAHs to derive star-formation rates in ULIRGs. The slopes of the observed correlations 
between PAH emission ratios and strengths and S'[9.7/im] are somewhat flatter than those 
expected if there is a screen of dust between the emission regions and the observer; that 
is, PAH emission is less absorbed than expected from the simplest model. This difference 
suggests that the PAH-emitting regions are spread out within the absorbing medium, rather 
than being concentrated in one central point. Because of the extinction, we are more likely 
to observe those of the PAH-emitting regions that are closest to us and are therefore the 
least obscured. 

On the contrary, as was previously concluded^, H2 is not affected by absorption. H2 
fluxes and ratios fail to show any correlation with the strength of silicate opacity (Figure 
E^-c), even though the S(3) line coincides in wavelength with the peak of the sihcate dust 
opacity curve (Figure [1]) and should be very strongly affected. In the absence of an observ- 
able dependence, the question arises just how strong a correlation would be expected if H2 
emission were behind the absorbing material. Molecular hydrogen is known to be present 
at a range of excitation temperatures^ and ortho-to-para ratios^, so even a simple model 
should take into account a wide range of physical conditions. To calculate the expected H2 
line ratios, I assume that dM oc T'^dT, where dM is the mass of H2 gas with excitation 
temperatures in the range T,T + dT. For each p, 1 calculate the emitted H2 spectrum for 
ortho-to-para ratios between 1 and 3. The range 2.5 < p < 5.0 is sufficient to reproduce 
the observed range of H2 excitation diagrams of galaxies (as shown in Supplementary Infor- 
mation). I then calculate H2 line ratios expected if the H2-emitting region were behind a 
screen of cold dust and select those of the line ratios that are expected to be most sensitive 
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to extinction. 

Even given a generous range of model conditions and intrinsic range of H2 line fluxes 
and ratios, their dependence on extinction would have been unambiguous if it were present 
in Figures [3]^ a-c). It is also noteworthy that models at zero extinction succeed in reproducing 
the observed range of line ratios, which confirms that the models adequately represent the 
range of physical conditions in real objects. Using different line combinations, more than a 
dozen line ratios can be computed; none shows dependence on S[9.7^in]. Neither the line 
ratios associated with the cooler component (median excitation temperature Texc = 330K 
as measured from S(l), S(2) and S(3)) nor the ones from the warmer component (median 
Texc = 1200K as measured from S(3), S(4) and S(7)) show any evidence for extinction. 

The fractional contribution of H2 and that of PAHs (extrapolated to low extinction) 
to the total luminosity is smaller in ULIRGs than in comparison galaxies. For PAHs, this 
finding is in agreement with the known non-linear relationship between bolometric and PAH 
luminosities^^. As for H2, there has to be a component directly associated with star forma- 
tion, just like in normal galaxies. However, as the extinction is increased, this component is 
quickly extinguished, and we can only see H2 emission that originates outside obscuration. 
The luminosity of this component may be estimated by observing that the H2/PAH ratios 
in ULIRGs exceed those seen in normal galaxies (Figure [3li). The data are well described by 
a model in which about as much (and up to 50% more) H2 emission is emitted outside the 
bulk of the obscuration as is directly associated with the starbursts, if the intrinstic PAH/H2 
ratios are the same in ULIRGs and in comparison galaxies. 

If dust is mixed in with emission regions or has patchy distribution, the correlations 
with silicate strength would be expected to be less steep than in the 'screen of dust' case. 
If dust is not completely cold, its thermal emission would fill in somewhat the absorption 
feature, making the observed 5'[9.7/im] an underestimate of the true optical depth, and the 
correlations would be expected to be steeper than in 'screen of dust' case. However, such 
effects cannot change the sign of the expected correlations. Since H2/PAH ratios increase 
with opacity when a decrease is expected from the models (Figure [3li), any model with the 
same spatial distribution of H2 and PAHs is ruled out by the data. 

The first step in understanding the origin of the rotational H2 emission in ULIRGs is the 
determination of the spatial extent of this emission, and in particular, whether this emission 
is confined to the outer parts of starburst regions or is extended on galaxy-wide scales. One 
way to proceed is to study the spatial distribution of the ro- vibrational lines in the near-IR, 
where observational capabilities for spatially resolved spectroscopy or narrow-band imaging 
are not as limited as in mid-IR. This approach yielded two successful observations of ULIRGs, 
one of NGC 6240^^ and one of Arp 220^^. In both cases the objects are composed of two 
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well-separated merging components, and the near-IR H2 emission peaks between them, near 
the colhsion front. 

Since most ULIRGs are either actively merging or dynamically disturbed^^, these obser- 
vations suggest the possibility that H2 lines trace shocks excited by inter-galactic interactions^^. 
Recent observations with Spitzer Space Telescope uncovered a class of unusually H2-luminous 
extragalactic objects^^'^^'^°, in which little or no star-formation is seen. One of these objects, 
a 40 kiloparsec inter-galactic shock in Stephan's Quintet^^'^^, emits tens of percent of its 
bolometric luminosity in rotational H2 lines. A similar process - albeit producing twenty 
times more mid-IR H2 luminosity in a median ULIRG than in Stephan's Quintet - may be 
operating in the outer parts of ULIRGs. Further support for this hypothesis is provided 
by observations of optical line emission from six ULIRGs^^, including four from this Letter, 
in which a contribution from inter-galactic shocks is directly seen. The median excitation 
temperature measured between the S(l) and S(3) transitions for ULIRGs with S[9.7fim]> 1 
in which the unobscured component dominates the observed H2 flux is 306 K, similar to that 
measured in Stephan's Quintet between the same transitions (350 K^^). 

Processes other than inter-galactic shocks could in principle lead to excitation of H2 
emission. However, as discussed in Supplementary Information, excitation by X-rays, by 
cosmic rays or by emission due to accretion onto the supermassive black hole are unlikely 
to be dominant producers of the observed H2 emission on energy grounds. Alternatively, H2 
could also be excited in shocks other than those produced by galactic collisions. For example, 
it could be due to super-galactic winds or supernova remnants which, having cleared up the 
environment around them, may be not as obscured as the photon-dominated regions that 
produce most of the PAH emission. These processes are neither supported nor ruled out 
as the dominant source of H2 excitation by the available data on ULIRGs. If they turn 
out to be important, this will be in contrast to lower- luminosity starburst galaxies in which 
these processes do not dominate H2 excitation^^. Future observations of spatial distribution 
and kinematics of H2 emission will help distinguish between intergalactic shocks and other 
scenarios. 

The observations of ULIRGs presented here raise the possibility that excitation of H2 in 
shocks and their subsequent cooling through the mid-IR rotational lines of H2 is a much more 
common phenomenon than previously thought. Instead of being an astronomical curiosity 
illustrated by a handful of H2-luminous objects and a few interacting galaxies, such cooling 
may be occurring in a wide range of conditions, as the large sample of ULIRGs demonstrates. 
If similar conditions are commonly encountered in the early universe^^, when H2 is one of the 
primary coolants, cooling via H2 emission in shocks can have a dramatic effect on formation 
of the first objects. 
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Supplementary Information is linked to tlie online version of the paper at 
www.nature.com/nature. 
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Supplementary Information 

ULIRGs have luminosities L > IO^^Lq, witli most of their emission coming out at 
infrared wavelengths^. The energy source of these objects remains a matter of controversy. 
Although the conventional interpretation is that they are powered by starbursts^, some 
authors argue that buried active galactic nuclei dominate energy output®. The IR spectra 
of ULIRGs display a wealth of features (Figure ISupp), including a dozen PAH emission 
features^'^", strong silicate absorption at lOyum and 18/im, an absorption complex at 6— 8/im 
due to water ice and hydrocarbons^^, as well as fine-structure atomic lines and rotational lines 
of H2. The measurements of features essential for the analysis of this Letter are described 
below and are listed in Tables 1 and 2 for the 48 ULIRGs from the 1 Jy sample®'^^ (hereafter 
the main sample). In addition, I use a more heterogeneous sample of 98 ULIRGs with silicate 
strength measurements^*^'^^ (hereafter the auxiliary sample), of which 68 objects have H2 
measurements in the literature'^. Twelve objects overlap between the main and the auxiliary 
samples. 

Among the comparison sample of SINGS galaxies^'^^'^°, I use nuclear pointings of 22 
star-forming galaxies, as well as 80 extra-nuclear pointings of star-forming galaxies and 
galaxies with active nuclei. In the latter case, the active nucleus is sufficiently far away 
that it does not affect the processes captured by the extra-nuclear pointings. Using these 
spectra, I obtain the ffux ratios of emission features shown with ellipses in Figures 2a,b 
and 3a,b,d. The distribution of these ratios is insensitive to whether only nuclear pointings 
or only extra-nuclear pointings or both are included. The H2/F(IR) ratios shown with 
ellipses in Figure 3c are obtained from nuclear pointings for star-forming galaxies, for which 
the IR ffuxes are calculated within the apertures of the spectroscopic observations^. The 
comparison PAH/F(IR) ratios shown with dotted ellipses in Figure 2c, d are obtained using 
stacked template spectra^*^ of all SINGS galaxies, without excluding those with active nuclei. 
Because of the uncertainties in the stacking procedure and because the set of galaxies included 
in this comparison is different from that used in other figures, the PAH/F(IR) ratios for 
galaxies should be regarded as an estimate. 

The apparent strength of silicate absorption is S[9.7fim] = — ln(/obs[9.7;um]//cont[9.7/im]), 
where fobs is the observed fiux density and /cont is the guess for the local continuum extrap- 
olated from pivot wavelengths outside the absorption feature. Defined this way, negative 
values of S[9.7^m] indicate silicate emission, while positive values indicate absorption, with 
S[9.7fira] ~ 4 for the most absorbed sources. To estimate /cont, I average emission in the 
wavelength ranges 5.3 — 5.6/im and 13.85 — 14.15/im and then interpolate between these 
points using a power-law. This method is almost identical to that used for the auxiliary 
sample^^; the minor difference at high optical depths is that I do not use a continuum 
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point at 7.7fim even in the cases of weak PAH emission. The strength of an absorption 
feature differs from the optical depth of this feature by an amount corresponding to the 
continuum opacity at this wavelength; for typical dust opacity curves in the literature, 
rdust[9.7/im] ~ (1.2 - 1.5) x S[9.7fim]. 

PAH[6.2/xm, 8.5/xm, ll.S/^m] fluxes are calculated by cutting out a < 3/im-wide part of 
the spectrum, excluding other features in that wavelength range, and then modeling this cut- 
out using a polynomial or a power-law continuum and Drude (damped harmonic oscillator) 
profiles with fixed widths^°. For PAH complexes the relative amplitudes of the components 
are fixed to their ratios in the template spectrum of normal star- forming galaxies^°. For 
example, within the 11.3/im complex the amplitude ratio of the 11.23/im and 11.33/im com- 
ponents is fixed to 1.25:1. Depending on the model for the local continuum, from a constant 
to a cubic polynomial, the number of fit parameters varies from two to five. Because the 
PAH[7.7/im] feature may be strongly affected by a number of absorption and emission fea- 
tures on both sides of the peak, it is fit over a narrow wavelength range (7.4— 8.0/im). Several 
measures of the luminosity of this feature are obtained by using different choices of the con- 
tinuum (constant, linear, and power-law), and the slope of the spectrum at 7.5/im is used 
as an additional luminosity measure after it is calibrated using the star-forming galaxies 
template^*^. The median value is used as a final flux. Example fits are shown in Figure 
HBupp. The quality of fits to PAH features is typically very good, and most are detected in 
almost all objects with high confidence. 

The strength of the 6/im absorption feature is defined similarly to the silicate strength 
and is computed by fitting a four-parameter model to the 5— 6/im cut-out from the spectrum, 

F^ = (lin. cont. + PAH[5.27;um] + PAH[5.70/im] + PAH[6.22/im])) exp {-S[6.0fim] x G(A)) , 

(1) 
where G{X) is an amplitude unity Gaussian centered at 6.0/im with a = 0.19/im. At 5/im < 

A < 6yum, G{X) is a good approximation to the laboratory opacity of water ice^^. At longer 

wavelengths, this fit is not appropriate, not only because G{X) no longer describes water ice 

opacity, but also because there may be other sources of opacity contributing to the 6— 8/xm 

complex^*^'^^ , plus the PAH feature at 6.2/im makes fitting this region very difficult. The 

PAH[5.27/im] and PAH[5.70/im] features in equation ([1]) are fit with two Drude profiles using 

a fixed 1:2 amplitude ratio and fixed widths^°. The PAH[6.22/xm] feature is included into the 

fit (with a fixed 15.7:1 ratio compared to PAH[5.27/im]), but it makes only a minor difference 

since only the weak extended wings of the profile contribute at wavelengths < 6/im. Water 

ice has little continuum opacity at these wavelengths, so that the measured strength of 

the feature is almost identical to the total optical depth, THjolS-O/xm] = 1.03 x ^[G.Oyum]. 

S'[6.0/im] is a rather difficult measurement because of the presence of multiple overlapping 

absorption and emission features in the same wavelength range. The relationship between 



silicate opacity and ice opacity is shown in Figure [SBupp. It appears that there may be an 
upper hmit to the amount of ice per given sihcate strength, ^[6.0/im] < 0.6 x S'[9.7/xm]. 
Deep sihcate absorption is a necessary but not sufficient condition for the presence of ice 
absorption. 

H2 hne fluxes are calculated by fitting Gaussian functions plus a quadratic continuum 
to cut-outs from the spectrum around the lines. The widths of the Gaussians are fixed to the 
resolution of the spectral orders that the lines fall into^°. The rotational transitions occur in 
the mid-IR and are denoted as S(J), where J is the angular momentum quantum number 
of the final state, into which the molecule transitions from the level J + 2. Observationally, 
S(l), S(2) and S(3) are the easiest to measure at low redshifts because of lack of blending 
with other features and relatively high signal-to-noise ratios within the spectral orders in 
question. Additionally, S(0), S(4) and S(7) can be measured in some cases. 

Observed line fluxes F{J) can be translated into level populations Nj: 

AttDIFj 

^J+2 - -. 7-^ ^TT, (2) 

where Dl is the luminosity distance to the source, A are the Einstein coefficients^'^ and 
Ejj^2 ~ Ej = hc/X is the energy of the transition. The rotational energy levels are given 

by34,35 

Ej = 85.35K • ksJi-J + 1) - 0.068K • keJ^iJ + 1)^ (3) 

where gj is the degeneracy of the levels and ks is the Boltzmann constant. For an equilibrium 
ortho-to-para (odd J-to-even J) ratio of 3, gj = 3(2 J + 1) for odd J and gj = 2J + 1 ioi 
even J. Any two or more line fluxes can be used to measure the excitation temperature as 
a ffiting parameter according to 

— (xexp\- . (4) 

gj V kBTe^cJ 

In practice, gas is present at a range of excitation temperatures, from ~150K to ~2000K. 
Typical two-temperature ffis to the extinction diagrams for ULIRGs^, such as those in Figure 
[EBupp, indicate that ~ IO^Mq of gas is found at T < 300K and ~ lO^M© of gas is found 
at T > lOOOK. The massive colder component contributes very little to lines with high J 
values, but dominates the fluxes of low- J lines, while a small amount of warmer gas boosts 
high- J fluxes. To quantify deviations from a single-temperature spectrum, for any three lines 
I can compute R[Ji/ J2, J3], which is the ratio of the observed flux of the S(Ji) line to that 
expected from lines S(J2) and 8(73) if all three transitions were from the same excitation 
temperature. In the case of multi-temperature gas at an equilibrium ortho-to-para ratio and 
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with no extinction, excitation diagrams are concave and R values are less than 1 (for Ji 
between J2 and J3). 

Bolometric fluxes are calculated by fitting IRAS or Spitzer broad-band data using a 
single-temperature black body function modified by emissivity oc u^^. The values obtained 
this way for IRAS data are within 10% of those obtained using the fitting formula 

F(IR) = 1.8 X 10-iierg/sec/cm'(13.48Fi2 + 5.I6F25 + 2.58^60 + Fioo), (5) 

where F12,... are IRAS fluxes in Jy^. 

The accuracy of PAH flux measurements is limited by the assumptions that were made 
about the shapes of the features which are fixed in my analysis. However, the results in 
this Letter appear to be robust to the fitting method used for calculating PAH fluxes. For 
example, using PAH fluxes obtained by summing up the spectrum above the local linear 
continuum® in Figures 2a, c produces similar results and the correlations with S'[9.7/im] are 
easily detectable. As for PAH[6.2/im] (Figure 2b), it needs to be corrected for water ice ab- 
sorption (i.e., multiplied by exp(S'[6.0;um]) before the correlation of PAH[11.3/im]/PAH[6.2yum] 
with S[9.7fim] becomes apparent. All conclusions remain the same if the main sample of 
ULIRGs is split by optical classification (LINERs vs HII regions) and each category is con- 
sidered separately. 

Previous studies of mid-IR H2 lines adopted varying views on the effects of absorption. 
For ULIRGs observed with ISO, H2 fluxes were corrected for absorption before constructing 
excitation diagrams^^. More recently^ zero absorption correction was suggested for H2 in 
the auxiliary sample on the basis of simultaneous measurements of S(l)/S(3) and S(l)/S(2) 
ratios as a function of the excitation temperature (this argument is analogous to saying that 
i?[2/l,3] values in ULIRGs agree with the model value). It is largely unknown if the near- 
IR ro-vibrational lines in ULIRGs are co-spatial with the mid-IR rotational lines; previous 
studies of ULIRGs in the near-IR tend to correct the entire spectrum for extinction^^. 

Several processes can lead to near-IR and mid-IR line emission from H2'^®. One possi- 
bility is excitation by UV photons of electronic levels in the molecule, which then radiatively 
decay into excited vibrational and rotational states. Other possibilities are collisional exci- 
tation in dense gas heated by UV or X-ray radiation, shocks or cosmic rays, or formation 
of hydrogen molecules in excited states on dust grains. Spitzer data by themselves can- 
not be used to establish the specific mechanism responsible for H2 excitation, because the 
low-energy transitions in question are almost certainly thermalized*^. 

Shock excitation is the explanation favored by this Letter; other possibilities appear 
unlikely. For example, ULIRGs presented here show no signs of an actively accreting super- 
massive black hole, in that no radio jets or high-ionization emission lines are seen in these 
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objects. Besides, any component of H2 emission excited by a deeply buried active nucleus 
would be heavily extincted. The observed X-ray fluxes are too smalF^ to be responsible 
for the observed H2 emission given a reasonable excitation efficiency of 10% for all H2 lines 
combined^*^. To estimate the contribution of cosmic rays, I assume a 10^" ergs of cosmic ray 
output from each supernova, and furthermore that 0.01 supernova is produced for every Mq 
of star formation and that the star-formation rate of ULIRGs is 100—500 M0/year^^. The 
luminosity in cosmic rays is then only a factor of 3—15 larger than the median luminosity of 
rotational H2 lines. Therefore, even if all cosmic rays are deposited directly into the inter- 
stellar medium of the galaxy (a scenario which is contradicted by some recent observations 
which suggest only 5% deposition efficiency'^^ ) , somehow avoiding regions with high density 
and thus high obscuration, tens of per cent of their energy would need to go directly into the 
observed H2 rotational emission. The combination of these requirements makes it unlikely 
that cosmic rays are the dominant source of H2 excitation. 

Table notes. S'[9.7/im] and 5'[6.0/xm] are the apparent strengths of the silicate and 
water absorption features, respectively. The small contribution of PAH and line emission 
to the wavelength ranges used in the calculation of /obs and /cont results in a systematic 
uncertainty in 5'[9.7/xm] which I take to be 0.21. This value is calculated from the variance 
of S'[9.7/im] measured in off- nuclear spectra of the comparison galaxies (when it is expected 
to be 0), while the statistical error is very small by comparison. I estimate the systematic 
uncertainty in S'[6.0/im] to be 0.16 by varying the functional form of the continuum (power 
law instead of linear) and by computing the variance of measured S'[6.0/im] for objects in 
which the S'[6.0/im] = fit is statistically acceptable. The nominal statistical uncertainty 
in the flt parameter is much smaller in most cases. Bolometric luminosities are calculated 
by fltting modifled black-body functions (see text). Uncertainty is taken to be 10%. PAH 
luminosities (L) and their uncertainties (a) are in Table 1. Measurements below 2(t are 
listed as 0.0. The uncertainties in PAH fluxes are calculated by adding in quadrature an 
estimate of the systematic error (obtained by varying the functional form of the continuum) 
and statistical error (derived from the standard deviation of the PAH amplitude in the linear 
flt). H2 line luminosities (L) and uncertainties {a) are in Table 2. For H2 fluxes, I use the 
standard deviation of the line amplitude in the linear flt to derive flux uncertainty. Fluxes 
are converted to luminosities using redshifts^ and an h = 0.7, Qm = 0.3, Q\ = 0.7 cosmology. 



REFERENCES 

1. Sanders, D.B. & Mirabel, I.F. Luminous infrared galaxies. Ann. Rev. Astron. Astrophys., 
34, 749-792 (1996) 



- 11 - 



2. Genzel, R. et al. What powers ultraluminous IRAS galaxies? Astrophys. J., 498, 579-605 

(1998) 

3. Higdon, S.J.U., Armus, L., Higdon, J.L., Soifer, B.T. & Spoon, H.W.W. A Spitzer Space 

Telescope Infrared Spectrograph survey of warm molecular hydrogen in ultraluminous 
infrared galaxies. Astrophys. J., 648, 323-339 (2006) 

4. Roussel, H. et al. Warm molecular hydrogen in the Spitzer SINGS galaxy sample. Astro- 

phys. J., 669, 959-981 (2007) 

5. Hollenbach, D.J. & Tielens, A.G.G.M. Dense photodissociation regions (PDRs). Ann. Rev. 

Astron. Astrophys., 35, 179-215 (1997) 

6. Appleton, P.N. et al. The dark side of reionization: probing cooling in the early universe. 

Astro2010: The Astronomy and Astrophysics Decadal Survey, Science White Papers, 



no. 2, |http://arxiv.org/abs/0903.1839| (2009) 



7. Houck, J.R. et al. The infrared spectrograph on the Spitzer Space Telescope. Astrophys. 



J. Supp., 154, 18-24 (2004). Data are available at |http://ssc. spitzer. caltech.edu 



8. Imanishi, M. et al. A Spitzer IRS low-resolution spectroscopic search for buried AGNs 

in nearby ultraluminous infrared galaxies: A constraint on geometry between energy 
sources and dust. Astrophys. J. Supp., 171, 72-100 (2007) 

9. Allamandola, L.J., Tielens, A.G.G.M. & Barker, J.R. Interstellar polycyclic aromatic 

hydrocarbons - The infrared emission bands, the excitation/emission mechanism, and 
the astrophysical implications. Astrophys. J. Supp., 71, 733-775 (1989) 

10. Hao, L. et al. The distribution of silicate strength of AGNs and ULIRGs. Astrophys. J., 

655, L77-L80 (2007) 

11. Spoon, H.W.W. et al. Mid-infrared galaxy classification based on silicate obscuration 

and PAH equivalent width. Astrophys. J., 654, L49-L52 (2007) 

12. Kennicutt, R.C., et al. SINGS: The SIRTF nearby galaxies survey. Proc. Astron. Soc. 



Pacific, 115, 928-952 (2003). Data available at http : / /sings . stsci . edu [ 



13. Shi, Y. et al. Aromatic features in AGNs: star-forming infrared luminosity function of 

AGN host galaxies. Astrophys. J., 669, 841 (2007) 

14. van der Werf, P.P. et al. Near-infrared line imaging of NGC 6240 - collision shock and 

nuclear starburst. Astrophys. J., 405, 522-537 (1993) 



-12- 



15. van der Werf, P.P. The nuclear interstellar medium of ultraluminous infrared galaxies. 

pp. 303-306 in Galaxy Interactions at Low and High Redshift, lAU proc. 186, eds. 
J.E.Barnes & D.B.Sanders (1999) 

16. Kim, D.-C, Veilleux, S. & Sanders, D.B. Optical and near-infrared imaging of the IRAS 

1 Jy sample of ultraluminous infrared galaxies. I. The atlas. Astrophys. J. Supp., 143, 
277-314 (2002) 

17. Rieke, G.H., Cutri, R.M., Black, J.H., Kailey, W.F., McAlary, C.W., Lebofsky, M.J., & 

Elston, R. IO^^Lq starbursts and shocked molecular hydrogen in the colliding galaxies 
Arp 220 (=IC 4553) and NGC 6240. Astrophys. J., 290, 116-124 (1985) 

18. Appleton, P.N. et al. Powerful high-velocity dispersion molecular hydrogen associated 

with an intergalactic shock wave in Stephan's Quintet. Astrophys. J., 639, L51-L54 
(2006) 

19. Egami, E., Rieke, G.H., Fadda, D. & Hines, D.C A large mass of H2 in the brightest 

cluster galaxy in Zwicky 3146. Astrophys. J., 652, L21-L24 (2006) 

20. Ogle, P. et al. Shocked molecular hydrogen in the 3C 326 radio galaxy system. Astrophys. 

J., 668, 699-707 (2007) 

21. Guillard, P., Boulanger, F., Pineau des Forets, G. & Appleton, P.N. H2 formation and 

excitation in the Stephan's Quintet galaxy-wide collision. Astron. Astrophys., 502, 
515-528 (2009) 

22. Monreal-Ibero, A., Arribas, S., & Colina, L. LINER-like extended nebulae in ULIRGs: 

shocks generated by merger- induced flows. Astrophys. J., 637, 138-146 (2006) 

23. Forbes, D.A., et al. High-resolution imaging of forbidden Fe II 1.64 microns, Brackett- 

gamma, and H2 1-0 S(l) emission in the starburst galaxy NGC 253. Astrophys. J., 
406, L11-L14 (1993) 

24. Greif, T.H., Johnson, J.L., Klessen, R.S., & Bromm, V. The first galaxies: assembly, 

cooling and the onset of turbulence. Mon. Not. R. Astron. Soc., 387, 1021-1036 
(2008) 

25. Weingartner, J.C. & Draine, B.T. Dust grain-size distributions and extinction in the 

Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud. Astrophys. J., 
548, 296-309 (2001) 

26. Kemper, F., Vriend, W.J. & Tielens, A.G.G.M. The absence of crystalline silicates in the 

diffuse interstellar medium. Astrophys. J., 609, 826-837 (2004) 



-13- 



27. Chiar, J.E. & Tielens, A.G.G.M. Pixie dust: the silicate features in the diffuse interstellar 

medium. Astrophys. J., 637, 774-785 (2006) 

28. Zakamska, N.L., Gomez, L., Strauss, M.A. & Krolik, J.H. Mid-infrared spectra of 

optically-selected type 2 quasars. Astron. J., 136, 1607-1622 (2008) 

29. Gerakines, P.A., Schutte, W.A., Greenberg, J.M. & van Dishoeck, E.F. The infrared 

band strengths of H2O, CO and CO2 in laboratory simulations of astrophysical ice 
mixtures. Astron. Astrophys.^ 296, 810-818 (1995) 

30. Smith, J.D.T. et al. The mid-infrared spectrum of star-forming galaxies: global properties 

of polycyclic aromatic hydrocarbon emission. Astrophys. J., 656, 770-791 (2007) 

31. Spoon, H.W.W. et al. Fire and ice: Spitzer Infrared Spectrograph (IRS) mid-infrared 

spectroscopy of IRAS F00183-7111. Astrophys. J. Supp., 154, 184-187 (2004) 

32. Kim, D.-C. & Sanders, D.B. The IRAS 1 Jy survey of ultraluminous infrared galaxies. 

I. The sample and luminosity function. Astrophys. J. Supp., 119, 41-58 (1998) 

33. Turner, J., Kirby-Docken, K., & Dalgarno, A. The quadrupole vibration-rotation transi- 

tion probabilities of molecular hydrogen. Astrophys. J. Supp., 35, 281-292 (1977) 

34. Huber, K.P. & Herzberg, G. Molecular spectra and molecular structure. IV. Constants 

of diatomic molecules. (Van Nostrand-Reinhold, New York, 1979) 

35. Landau, L.D. & Lifshitz, E.M. Quantum mechanics: non-relativistic theory. 

(Butterworth-Heinemann - Oxford, Boston, 1991) 

36. Rigopoulou, D., Kunze, D., Lutz, D., Genzel, R., & Moorwood, A.F.M. An ISO-SWS 

survey of molecular hydrogen in starburst and Seyfert galaxies. Astron. Astrophys., 
389, 374-386 (2002) 

37. Goldader, J.D., Joseph, R.D., Doyon, E., & Sanders, D.B. Spectroscopy of luminous 

infrared galaxies at 2 microns. II. Data for galaxies with 11.2 < log(LiR/L0) < 11.9. 
Astrophys. J. Supp., 108, 449-470 (1997) 

38. Shull, J.M. & Beckwith, S. Interstellar molecular hydrogen. Ann. Rev. Astron. Astrophys., 

20, 163-190 (1982) 

39. Teng, S.H., Wilson, A.S., Veilleux, S., Young, A.J., Sanders, D.B., & Nagar, N.M. A 

Chandra X-ray survey of ultraluminous infrared galaxies. Astrophys. J., 633, 664-679 
(2005) 



-14- 



40. Lepp, S. & McCray, R. X-ray sources in molecular clouds. Astrophys. J., 269, 560-567 

(1983) 

41. Dopita, M.A. et al. Modeling the pan-spectral energy distribution of starburst galaxies. 

Astrophys. J., 619, 755-778 (2005) 

42. Acero, F., et al. Detection of gamma rays from a starburst galaxy. Science., 326, 1080- 

1082, (2009) 

43. Armus, L. et al. Observations of ultraluminous infrared galaxies with the Infrared Spec- 

trograph (IRS) on the Spitzer Space Telescope: early results on Markarian 1014, 
Markarian 463, and UGC 5101. Astrophys. .J. Supp., 154, 178-183 (2004) 

44. Spoon, H.W.W. et al. Ice features in the mid-OR spectra of galactic nuclei. Astron. 

Astrophys., 385, 1022-1041 (2002) 

45. Armus, L. et al. Detection of the buried active galactic nucleus in NGC 6240 with the 

Infrared Spectrograph on the Spitzer Space Telescope. Astrophys. J., 640, 204-210 
(2006) 



This preprint was prepared with the AAS E^TJtjX macros v5.0. 



-15- 




10 , ,14 

wavelength (fim) 



Fig. 1. — Wavelengths of emission features present in ULIRG spectra and repre- 
sentative opacity curves. If emission regions are embedded in dust or ice, those emission 
features that are near the peaks of opacity should be the ones most strongly affected. If H2 
emission originates inside silicate dust obscuration (example opacities^^'^^'^^ in black), then 
S(3) should be strongly extincted, S(l) somewhat less so, and other lines less still. If PAH 
emission originates inside obscuration, the features centered at 8.5/im and 11.3/im should be 
more affected by silicates than the other ones, while PAH[6.2/im] is the only feature that may 
be affected by water ice (grey; smoothed data for the whole complex^^ and laboratory data 
for water ice^^, shown at the maximum strength according to the 5'[6.0/xm]= 0.6S'[9.7;Um] 
relation - see Supplementary Information). 
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Fig. 2. — PAH features are affected by dust extinction. This conclusion is evidenced 
by correlations of the observed ratios of PAH fluxes (a,b) and the PAH/F(IR) ratios (c,d) 
with the apparent strength of the silicate opacity feature. Spearman rank probabilities 
of the null hypothesis that the plotted values are uncorrelated are (a) P[NH]= 10"*^, (c) 
P[NH]= 10-^ (d) P[NH]= 10*^ In panel (b), the observed PAH[11.3/im]/PAH[6.2/im] 
ratio (grey circles) is uncorrelated with the absorption strength (P[NH]= 0.22), but when 
PAH[6.2/im] is corrected for water ice absorption (black circles), the correlation becomes 
apparent (P[NH]= 2 x 10^^), suggesting that the PAH-emitting region is located behind both 
silicate absorption and water ice absorption. In each panel (a) and (b), PAH ratios calculated 
for the comparison sample of nearby star-forming galaxies^^'^'^ are shown by an ellipse whose 
semi-axes are determined by the standard deviation of the corresponding measure. Although 
PAH ratios are known to vary as a function of physical conditions'^, PAH ratios seen in 
ULIRGs are consistent with those found in the comparison sample when extrapolated to 
low absorption. Dashed and dotted lines illustrate how unobscured ratios would change in 
the presence of an increasing amount of cold dust between the emitter and the observer for 
representative opacity curves^^'^^. In panels (c) and (d), the model calculation assumes that 
all absorbed flux is re-emitted at longer wavelengths, but the total flux does not change. 
PAHs constitute a higher fraction of the total luminosity output in low-luminosity galaxies 
(estimates shown with dotted ellipses) than they do in ULIRGs^'. The grey shaded areas 
show the la range in the vertical offset around the best linear fit for ULIRG data. Filled 
circles denote detections, arrows denote 3a upper limits, and la standard errors are shown 
in each panel except (b) where individual error bars are omitted for clarity and the median 
error is shown in the bottom left corner. 
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Fig. 3. — H2 emission in ULIRGs is not affected by extinction and shows excess 
over the H2/PAH ratio observed in normal galaxies. The ratios of H2 fluxes (a,b) and 
the H2/F(IR) ratios (c) for ULIRGs from two samples (circles^, crosses^) are uncorrelated 
with the apparent strength of the silicate absorption. In (b), R[Ji/ J2, J3] is the ratio of 
the observed flux of the line S(Ji) to that expected based on S(J2) and 8(73) assuming all 
three lines come from a region with a single excitation temperature. Grey areas in (a, b) 
and lines in (c) show the expected trends of line ratios with apparent silicate strength if H2 
emission is behind a screen of dust. Models in dark grey assume an ortho-to-para ratio of 3 
and a realistic range of excitation temperatures, whereas light grey areas are an extension 
of the model to include ortho-to-para ratios between 1 and 3. Although correlations are 
expected if H2 is affected by silicate absorption, none are detected (P[NH]= 5 — 75%). In 
(d), if H2 and PAH emission had the same spatial distribution, the H2/PAH ratios would 
be expected to decrease (dotted and dashed lines) because dust opacity at the wavelength 
of the H2 line is greater than that at the wavelelngth of the PAH feature, but in fact an 
increase is observed (P[NH] =0.007, 0.016 for top and bottom). The H2/PAH ratios are 
described better by a model that combines an obscured component of H2 associated with 
star formation and an unobscured H2 component with the luminosity equal or somewhat 
greater than that of the buried component (solid lines in panels (d), top and bottom as- 
sume L[H2, outer] /L[H2, inner] = 1 and 1.5, respectively). Arrows denote 3a upper and 
lower limits, and la standard errors are shown for all measurements. Ellipses show flux 
ratios for comparison galaxies (semi-axes are determined by the standard deviation of the 
corresponding measure) . 



-18- 



c 

T3 
X 

2 0.1 



O.B - 



Arp 820 (IRAS 15327 + 2340) 




10 15 

rest — frame wavelength, /j.m 



Fig. 4. — Supp. Example fits of emission and absorption features. The black dotted 
line shows the power-law continuum used for determination of the strength of the silicate 
feature (black arrow). Fit to water ice absorption is shown in grey, with the dotted line for 
the continuum and the arrow for the strength determination. Drude fits to PAH emission 
features are in light blue and Gaussian fits to H2 emission lines are in dark blue. 
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Fig. 5. — Supp. Relationship between the strengths of dust and water ice absorption 
features. Triangles show measurements for the main sample and stars show UGC 5101^^, 
IRAS 00183-7111^^, and objects from the auxiliary sample with 5'[6.0//m] measurements'^^. 
Among the latter, only the four sources with wavelength coverage starting at A < 5.6/xm were 
included (in the absence of short-wavelength coverage it is difficult to distinguish between ice 
absorption and PAH emission). The dotted line shows S'[6.0/.im] = 0.6S[9.7fini\ and the error 
bars in the upper right corner represent estimated systematic uncertainties of the feature 
strength measures. 
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Fig. 6. — Supp. Example H2 excitation diagrams (black) for objects with five or more 
rotational lines detected, including several ULIRGs from the main sample marked 'IRAS' 
and several star-forming galaxies'^'^^ marked 'NGC. Objects are sorted by the S(l)/S(7) ratio, 
with lower ratios on top. The observer assumes an equilibrium ortho-to-para ratio of 3 to 
compute values on the vertical axis from line fluxes (eq. [2]). Excitation temperature is defined 
as the inverse of the slope of the curves, so high- J transitions tend to have higher excitation 
temperatures. Arrows mark energies of the upper level with Jupper which radiatively decay 
into the level J = Jupper — 2; the corresponding line is denoted S{J). Thick grey curves show 
excitation diagrams for models discussed in the Letter (gas emits at a range of excitation 
temperatures, with dM oc T~^dT). The parameter p increases from 2.5 to 5.0 from top to 
bottom. Lighter zigzag lines show model excitation diagrams for gas with intrinsic ortho-to- 
para ratio of 1. Some objects show deviations from a concave excitation diagram expected in 
the case of the equilibrium ortho-to-para ratio. For each individual object, it is not possible 
to distinguish between a non-equilibrium ortho-to-para ratio (such as those seen in some 
nearby galaxies^) and extinction on the basis of excitation diagrams alone, since the same 
lines (S(l) and S(3)) may suffer from both effects. The lack of a correlation between line 
ratios and S'[9.7/im] or line deficits (parametrized by R[Ji/J2, J3]) and S[9.7fim\ either for 
the whole sample of ULIRGs or for the well-measured subsample shown here suggests that 
absorption is not important. 



Table 1. Absorption strengths and PAH luminosities in ULIRGs 
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S[9.7Mm] 


S[6.0/Jm] 


L[IR] 


L[6.2Mm] 


cr[6.2^m] 


L[7.7fim] 


a[7.7^ni] 


i[8.5Mml 


cr [8. 5 fiin] 


L[11.3Mm] 


(T[11.3/^m] 








10*-'' crg/s 


10''^ crg/s 


10*^ org/s 


10*^ crg/s 


Uf^ crg/s 


10*^ crg/s 


10*^ crg/s 


10"^ crg/s 


10*^ crg/s 


IRAS F00091 - 0738 


2.86 


0.83 


7.64 


6.9 


1.4 


47.1 


5.4 


0.0 


1.8 


5.2 


0.2 


IRAS F00188 - 0856 


2.60 


1.14 


9.66 


10.7 


1.0 


69.4 


4.3 


15.5 


2.0 


7.8 


0.3 


IRAS F00456 - 2904 


0.74 


0.26 


6.77 


26.9 


1.6 


98.1 


3.9 


29.7 


1.9 


16.4 


0.4 


IRAS F00482 - 2721 


1.84 


0.51 


4.50 


4.9 


0.8 


14.7 


3.9 


0.0 


2.0 


2.5 


0.2 


IRAS F01004 - 2237 


0.44 


0.12 


7.71 


6.6 


1.4 


69.0 


4.7 


0.0 


2.3 


3.5 


0.4 


IRAS F01166 - 0844 


2.73 


0.23 


4.98 


2.2 


1.1 


22.8 


3.8 


0.0 


1.7 


2.8 


0.3 


IRAS F01298 - 0744 


3.43 


0.89 


9.50 


4.2 


1.4 


88.1 


6.7 


0.0 


3.3 


5.1 


0.4 


IRAS F01569 - 2939 


2.17 


0.48 


7.25 


6.6 


2.1 


43.4 


5.4 


0.0 


3.2 


6.5 


0.5 


IRAS F02411 + 0353 


0.62 


0.00 


6.71 


34.7 


2.5 


142.5 


6.3 


35.1 


2.9 


19.5 


0.5 


IRAS F03250 + 1606 


1.15 


0.00 


5.14 


18.8 


1.2 


64.2 


4.7 


17.8 


2.1 


10.1 


0.3 


IRAS F04103 - 2838 


0.51 


0.18 


6.28 


20.3 


1.2 


83.6 


5.4 


22.2 


2.5 


14.0 


0.3 


IRAS F08572 +3915 


3.58 


0.12 


5.12 


0.0 


1.4 


0.0 


3.6 


0.0 


1.7 


0.0 


0.4 


IRAS F09039 + 0503 


1.83 


0.96 


5.20 


7.2 


0.9 


28.2 


4.1 


0.0 


1.9 


4.9 


0.3 


IRAS F09116 + 0334 


0.76 


0.00 


5.76 


22.5 


1.6 


98.5 


3.9 


28.0 


2.6 


16.3 


0.5 


IRAS F09539 + 0857 


3.03 


0.56 


3.77 


6.4 


1.0 


31.7 


3.9 


0.0 


2.1 


3.7 


0.2 


IRAS F10190 + 1322 


0.76 


0.00 


4.49 


21.2 


1.4 


81.0 


3.8 


24.9 


1.1 


14.0 


0.3 


IRAS F10378 + 1108 


2.41 


1.30 


8.88 


4.2 


1.3 


47.0 


4.0 


0.0 


2.4 


4.5 


0.3 


IRAS F 10485 - 1447 


2.67 


0.84 


6.79 


8.3 


1.3 


38.8 


3.8 


7.3 


2.1 


3.9 


0.2 


IRAS F10494 + 4424 


1.24 


0.36 


6.61 


15.7 


1.2 


65.4 


4.1 


16.5 


1.3 


6.9 


0.2 


IRAS F11095 - 0238 


3.13 


0.64 


7.50 


4.4 


1.0 


42.4 


4.7 


0.0 


1.7 


4.4 


0.3 


IRAS F11130 - 2659 


2.65 


0.59 


5.10 


5.3 


1.4 


38.0 


3.9 


0.0 


2.3 


3.4 


0.2 


IRAS F11387 + 4116 


0.69 


0.00 


5.33 


14.1 


1.3 


51.9 


4.3 


14.7 


2.7 


9.3 


0.3 


IRAS F11506 + 1331 


2.06 


0.50 


9.19 


21.8 


1.5 


86.2 


4.5 


19.9 


2.0 


8.0 


0.3 


IRAS F12112 +0305 


1.34 


0.88 


9.09 


14.1 


1.1 


55.9 


4.5 


12.2 


1.0 


7.2 


0.3 


IRAS F12127 - 1412 


2.29 


0.37 


5.87 


3.3 


1.0 


38.7 


4.8 


0.0 


2.1 


3.0 


0.3 


IRAS F12359 - 0725 


1.97 


0.63 


5.43 


7.3 


1.2 


31.6 


4.9 


6.6 


2.3 


3.4 


0.2 


IRAS F13335 - 2612 


0.93 


0.24 


5.02 


17.2 


0.9 


64.5 


3.8 


17.4 


2.0 


11.0 


0.3 


IRAS F13509 + 0442 


1.14 


0.00 


6.95 


18.6 


2.1 


81.4 


6.6 


23.1 


2.3 


10.3 


0.3 


IRAS F13539 + 2920 


1.09 


0.27 


4.79 


19.5 


1.1 


75.8 


4.1 


18.2 


1.7 


9.4 


0.3 


IRAS F14060 + 2919 


0.60 


0.00 


5.03 


23.3 


1.7 


98.9 


4.2 


27.1 


2.4 


13.9 


0.4 


IRAS F14252 - 1550 


0.88 


0.00 


6.24 


9.9 


1.4 


43.9 


4.0 


10.3 


2.9 


5.5 


0.4 


IRAS F 14348 - 1447 


1.86 


1.12 


9.46 


14.6 


1.2 


52.7 


4.2 


8.5 


1.0 


8.2 


0.3 


IRAS F15206 + 3342 


0.47 


0.00 


6.41 


29.7 


1.9 


109.8 


5.5 


26.4 


3.0 


17.4 


0.5 


IRAS F15225 + 2350 


2.02 


0.69 


5.70 


12.3 


1.6 


57.6 


4.6 


10.0 


2.4 


6.6 


0.3 


IRAS F15327 + 2340 


2.69 


1.11 


6.17 


3.7 


0.7 


19.2 


3.6 


1.9 


0.3 


2.0 


0.2 


IRAS F16090 - 0139 


2.16 


0.85 


14.59 


19.7 


1.7 


103.1 


5.9 


0.0 


3.1 


14.4 


0.5 


IRAS F16468 + 5200 


2.04 


0.84 


5.04 


10.6 


1.5 


35.8 


4.9 


0.0 


3.0 


4.1 


0.3 


IRAS F16474 + 3430 


1.08 


0.24 


6.06 


22.8 


1.3 


100.4 


4.9 


25.0 


2.2 


12.6 


0.3 


IRAS F16487 + 5447 


1.41 


0.78 


6.44 


7.5 


1.1 


32.2 


4.4 


6.9 


1.5 


5.2 


0.2 


IRAS F17028 + 5817 


1.15 


0.53 


6.30 


10.5 


1.1 


51.9 


3.9 


14.3 


1.5 


5.8 


0.2 


IRAS F17044 + 6720 


1.62 


0.14 


5.70 


7.0 


1.2 


46.6 


5.2 


9.4 


2.4 


4.2 


0.4 


IRAS F20414 - 1651 


1.57 


0.53 


6.68 


7.7 


0.9 


33.5 


3.7 


8.1 


0.9 


3.9 


0.2 


IRAS F21208 - 0519 


0.96 


0.00 


4.49 


11.8 


1.2 


49.3 


3.9 


13.9 


2.0 


6.2 


0.2 


IRAS F21329 - 2346 


1.94 


0.66 


5.73 


8.8 


0.9 


36.0 


4.4 


5.7 


1.8 


3.5 


0.2 


IRAS F22206 - 2715 


1.28 


0.88 


6.86 


10.9 


1.3 


40.0 


4.9 


9.1 


2.3 


6.0 


0.3 


IRAS F22491 - 1808 


1.12 


0.00 


6.12 


10.3 


1.0 


40.6 


3.7 


9.4 


0.9 


6.2 


0.2 


IRAS F23234 + 0946 


1.77 


0.49 


5.68 


8.3 


1.0 


35.3 


3.9 


7.7 


1.9 


4.7 


0.2 


IRAS P23327 + 2913 


1.08 


0.00 


5.27 


4.6 


1.0 


20.8 


3.6 


6.3 


1.4 


3.4 


0.2 



to 

o 



Table 2. Ho luminosities in ULIRGs 



Object ID 


L[S(1)] 


^[S(l)] 


L[S(2)1 


<t[S(2)] 


L[S(3)1 


<t[S(3)] 


L[S(4)] 


<t[S(4)1 


Z,[S(7)1 


<t[S(7)] 




lO*"" org/s 


10*" org/s 


10*" crg/s 


10*" org/s 


10*" org/s 


10*" crg/s 


10*" org/s 


10*" org/s 


10*" crg/s 


10*" org/s 


IRAS F00091 - 0738 


47.2 


3.6 


0.0 


5.2 


26.3 


2.9 


0.0 


12.9 


17.5 


6.2 


IRAS F00188 - 0856 


76.2 


11.4 


17.2 


6.8 


23.5 


4.3 


0.0 


17.3 


0.0 


16.3 


IRAS F00456 - 2904 


33.4 


2.3 


20.1 


3.0 


39.3 


3.9 


0.0 


19.3 


25.3 


10.5 


IRAS F00482 - 2721 


0.0 


2.9 


0.0 


4.5 


0.0 


4.5 


0.0 


9.4 


0.0 


6.0 


IRAS F01004 - 2237 


47.6 


11.7 


36.7 


12.0 


48.6 


8.9 


0.0 


30.1 


0.0 


9.8 


IRAS F01166 - 0844 


23.7 


5.2 


11.8 


3.6 


0.0 


3.6 


11.9 


5.3 


0.0 


9.5 


IRAS F01298 - 0744 


60.7 


16.5 


13.3 


4.4 


13.9 


2.6 


48.7 


19.0 


0.0 


11.3 


IRAS F01569 - 2939 


139.9 


8.8 


28.0 


8.7 


83.8 


10.8 


0.0 


14.6 


0.0 


29.8 


IRAS F02411 + 0353 


116.3 


19.6 


0.0 


7.4 


33.6 


7.0 


0.0 


29.5 


0.0 


26.9 


IRAS F03250 + 1606 


0.0 


7.0 


12.9 


1.1 


24.1 


6.6 


0.0 


20.5 


0.0 


14.7 


IRAS F04103 - 2838 


63.2 


10.9 


0.0 


8.7 


11.1 


1.6 


0.0 


21.1 


33.8 


9.5 


IRAS F08572 + 3915 


0.0 


7.3 


24.5 


6.9 


9.5 


2.7 


0.0 


19.0 


0.0 


11.0 


IRAS F09039 + 0503 


149.2 


10.0 


48.5 


2.8 


76.1 


3.1 


21.8 


8.6 


44.8 


6.9 


IRAS F09116 + 0334 


0.0 


15.8 


0.0 


22.2 


11.9 


2.8 


0.0 


21.1 


0.0 


20.3 


IRAS F09539 + 0857 


55.9 


11.4 


10.8 


3.1 


20.2 


3.4 


16.2 


4.3 


28.0 


10.1 


IRAS F10190 + 1322 


52.8 


3.6 


7.9 


2.1 


21.1 


3.6 


0.0 


23.0 


22.2 


6.3 


IRAS F10378 + 1108 


57.3 


17.8 


23.5 


4.7 


42.7 


6.4 


26.9 


8.0 


158.9 


23.9 


IRAS F10485 - 1447 


30.4 


11.7 


14.2 


5.6 


13.6 


3.9 


0.0 


10.2 


38.2 


10.8 


IRAS F1D494 + 4424 


44.3 


3.3 


10.2 


3.6 


16.6 


3.2 


0.0 


18.1 


0.0 


10.1 


IRAS F11095 - 0238 


95.2 


6.9 


29.1 


1.7 


48.2 


2.7 


0.0 


15.1 


36.3 


5.5 


IRAS F11130 - 2659 


35.9 


6.8 


0.0 


4.9 


13.2 


4.5 


0.0 


7.4 


0.0 


17.8 


IRAS F11387 + 4116 


65.6 


11.5 


14.2 


3.6 


0.0 


4.4 


0.0 


18.7 


0.0 


15.2 


IRAS F11506 + 1331 


22.8 


8.9 


19.1 


5.1 


22.4 


2.6 


0.0 


25.7 


18.6 


6.8 


IRAS F12112 + 0305 


72.2 


10.0 


20.4 


3.4 


34.0 


4.7 


0.0 


19.2 


0.0 


23.1 


IRAS F12127 - 1412 


0.0 


10.6 


16.4 


3.4 


0.0 


7.4 


0.0 


8.2 


0.0 


17.3 


IRAS F12359 - 0725 


52.0 


13.8 


0.0 


8.4 


12.3 


5.1 


0.0 


9.6 


0.0 


16.2 


IRAS F13335 - 2612 


110.6 


8.1 


22.8 


2.7 


46.9 


3.8 


0.0 


19.0 


0.0 


11.4 


IRAS F13509 + 0442 


46.6 


11.9 


12.3 


5.0 


38.7 


9.4 


0.0 


22.0 


0.0 


22.9 


IRAS F13539 + 2920 


120.9 


3.5 


25.6 


4.9 


42.6 


2.4 


0.0 


21.3 


0.0 


12.2 


IRAS F14060 + 2919 


43.8 


7.7 


0.0 


5.9 


24.3 


4.2 


0.0 


30.6 


0.0 


21.1 


IRAS F14252 - 1550 


0.0 


10.5 


0.0 


12.0 


26.9 


6.0 


0.0 


14.4 


0.0 


17.3 


IRAS F 14348 - 1447 


109.3 


21.1 


32.9 


3.0 


50.2 


4.2 


0.0 


20.2 


52.3 


14.2 


IRAS F15206 + 3342 


86.3 


12.3 


0.0 


7.4 


34.6 


9.2 


0.0 


38.2 


0.0 


17.7 


IRAS F15225 + 2350 


0.0 


11.3 


0.0 


12.3 


21.0 


5.6 


0.0 


21.2 


0.0 


13.7 


IRAS F15327 + 2340 


19.7 


3.3 


5.7 


0.3 


7.9 


0.3 


0.0 


4.3 


19.0 


2.3 


IRAS F16090 - 0139 


122.8 


12.1 


42.1 


6.9 


65.0 


4.8 


82.2 


19.6 


109.4 


19.6 


IRAS F16468 + 5200 


30.6 


5.9 


0.0 


6.2 


19.7 


5.8 


0.0 


17.3 


0.0 


22.2 


IRAS F16474 + 3430 


72.6 


4.1 


18.2 


5.2 


22.4 


4.6 


0.0 


36.3 


0.0 


14.1 


IRAS F16487 + 5447 


95.8 


5.1 


16.1 


3.3 


58.7 


7.7 


19.9 


9.7 


0.0 


11.7 


IRAS F17028 + 5817 


46.7 


6.3 


5.2 


1.5 


11.9 


3.2 


0.0 


20.0 


0.0 


11.2 


IRAS F17044 + 6720 


52.5 


2.9 


0.0 


14.0 


20.6 


8.0 


0.0 


15.3 


0.0 


23.0 


IRAS F20414 - 1651 


39.2 


9.1 


0.0 


3.0 


6.3 


1.6 


0.0 


7.7 


0.0 


8.5 


IRAS F21208 - 0519 


39.4 


7.2 


9.1 


3.5 


13.3 


5.5 


0.0 


14.1 


0.0 


10.6 


IRAS F21329 - 2346 


0.0 


15.0 


11.1 


4.0 


31.0 


3.0 


0.0 


14.3 


0.0 


15.9 


IRAS F22206 - 2715 


55.5 


5.9 


0.0 


12.9 


36.9 


9.2 


0.0 


7.5 


0.0 


15.2 


IRAS F22491 - 1808 


44.3 


6.5 


9.0 


3.0 


15.4 


2.5 


0.0 


12.2 


32.2 


5.9 


IRAS F23234 + 0946 


98.2 


13.6 


19.0 


3.7 


49.3 


6.8 


0.0 


11.8 


0.0 


13.8 


IRAS P23327 + 2913 


58.3 


5.9 


12.7 


4.3 


19.3 


2.4 


0.0 


5.7 


0.0 


10.1 



